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1.  INTRODUCTION 


In  life  testing  or  survival  studies,  the  observation  of  the  time  of 
occurrence  of  a  failure  or  death  uv  be  prevented  bv  the  occurrence  of 
•  >>se  other  event,  resulting  in  a  loss  of  an  Item  or  Individual  from  the 
study.  In  this  type  of  censoring,  only  the  tins*  of  loss  can  be  observed 
when  the  loss  occurs  before  the  death  of  the  Item.  The  problem  of 
nonparaaetrlcal 1 v  estivating  a  survival  function  Iron  such  censored  data 
has  received  auch  attention  in  the  recent  statistical  literature.  Rreslow 
4  Crowley  (1974)  and  l.agakos  (1979)  have  given  excellent  reviews  of  this 


subject  . 

©  •  © 

Spec  If  leal  1)  ,  wc  consider  the  following.  Let  X,,  X, . X  be  the 

l  i  n 


true  survival  tines  of  n  it  res  or  individuals  which  are  censored  frc>«s  the 


right  by  a  sequen. e  L  .  V, . I  which  any  be  either  constants  or  random 

I  *  n 


variables.  It  is  aaauaed  chat  the  are  Independent,  Identically 


distributed  random  variables  with  a  common  unknown  distribution  function 
F(t).  We  wish  to  estimate  the  survival  function  F(t)  "  1  -  F(t)  •  Pr{X  >t) 
baaed  on  observations  consisting  of  a  sequent e  of  pairs  (Xj.dj),  where 


j  1  if  Xj  S 

Xi  “  oin  Xi,lr  an<)  ii  *  1  0  If  X*  >  Ult  1  •  1,  2, 


i 


1  I 


,  n. 


Thus,  ft  Is  known  which  observations  are  times  of  deaths  and  which  ones 
ar»  times  of  losses  (censored). 

One  of  the  most  popular  estimators  of  T(t)  Is  the  product  limit 
estimator  P(t)  proposed  by  Kaplan  4  Meier  (191X8),  which  was  shown  to  be 
self-consistent  bv  Efron  (1967).  Rreslow  4  Crowley  (1974)  gave  a  rigorous 
derivation  of  the  large  sample  properties  of  the  product  limit  estimator. 


Another  method  of  deriving  the  product  limit  estimate  In  a  maximum 
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likellhood  framework  considered  by  Nelson  (1969)  end  Brealow  (1972,  1974) 
was  to  restrict  the  set  of  distribution  functions  to  those  having  a  hazard 
function  which  was  constant  between  the  distinct  uncensored  failure  timet. 

The  product  Unit  estimate  Is  a  step  function  and  Is  not  well-defined  when 
the  largest  observation  Is  a  loss.  To  improve  on  this  situation,  Susarla 
4  Van  Ryzlr  (1976,  1978)  and  Ferguson  4  Ph.*dla  (1979)  proposed  nonparametrlc 
Bayes  estimators  which  are  defined  everywhere,  use  all  of  the  censored  and 
uncensored  observations,  and  result  in  smoother  estimates  than  the  product 
limit  estimate.  Susarla  &  Van  Rvzln's  (1976)  estimator  still  has  jumps 
at  the  uncenaored  (allure  times,  however. 

Cox  (1972,  1975)  proposed  the  proportional  hazard  model  with  covarlate 
variables  and  used  the  partial  likelihood  principle  to  analyze  survival 
data . 

lr.  many  situations,  the  life  distribution  Fit)  mav  be  assumed  or  known 
to  have  a  monotone  hazard  or  failure  rate  function  (Barlow  4  Proschan  (1975)). 
For  uncensored  data ,  Cr cnander  (1956)  derived  mexlmum  livelihood  estimators 
for  the  failure  rate  function  r(t)  and  F(t)  assuming  only  that  r(t)  Is 
Increasing.  (Throughout  the  paper  we  write  "increasing”  for  "nondecreaelng" . ) 
The  consistency  and  asymptotic  distribution  of  this  maximum  likelihood 
estimator  were  established  bv  Marsftall  4  Proschan  (1965)  anJ  Prakasa  Rao 
(190),  respectively. 

In  thle  paper  ve  obtain  the  maximum  likelihood  estimator  F(t)  of  F(t) 
under  th«  condition  t ha t  Fit)  has  a  monotone  failure  rate  baaed  on  a  set 
of  observations  subject  to  arbitrary  right  censorship  as  described 
previously.  This  maximum  likelihood  estljritor  Is  continuous  and  well-defined 
for  all  t  ?  0.  In  Section  3  an  example  Is  presented,  end  In  Section  4  en 


Indication  of  the  small  sample  properties  of  F(t)  as  compared  with  P(t)  Is 
given,  resulting,  f roc  a  Monte  Carlo  study  for  several  Welbull  distributions 
with  Increasing  failure  rate  functions.  The  estlaator  F(t)  performs  very 
well  in  the  tails  of  the  distribution  and  for  samples  under  rather  severe 
censorship. 


2.  THE  MAXIMUM  L I KX1. 1  HCX)l)  ESTIMATOR 

Let  1-1 . n  denote  the  sample  described  In  Section  1, 

and  let  f  denote  the  common  probability  density  function  of  the  X^.  Assume 

that  l',  ....  U  are  either  constants  or  Independent  randoe.  variables 
1  n 

which  are  also  Independent  of  Xj .  ....  X^.  Then  the  likelihood  function 
can  be  written  as  (l-a&akos  (1979)) 

n  d  _  1-d 

L  •  L\(x  ,d  ):  1  -  I,  ....  n?  •  "  (f(x  ))  'F(x  )) 

1-1 

Let  the  failure  rate  function  be  r(t)  -  f(t)/F(t).  It  follows  that  the 
likelihood  function  Is 

n  61  - 

L  ■  I!  (r(x  ) 1  Fix  ),  (2.1) 

1-1 

and  since 

F(t)  -  exp ,-/q  r(u)du),  t  t  0,  (2.2) 

we  can  write  from  (2.1) 

n  n  x 

In  L  •  I  d  In  r(x.)  ~  I  /n  r(u)du.  (2.3) 

1-1  1  1-1  0 

Without  any  loea  of  generality,  assume  that  Xj  &  ■  ...  S  x^.  We  also 


aasune  that  r(t)  Is  increasing.  Consider  the  failure  rate  function 


0. 


r*(t)  - 


rUj).  xj  '  *  x1+1»  1 

r  (x  ) .  *  s  i  • 
n  n 


...  n-1 


Then  for  each  t.  r(t)  2  r*U),  ami  from  (2.3),  we  obtain 


In  L  s  I  i.  In  r(x  )  -  1  f0  r*iu)du 

1-1  1-1 


-  £  In  rl*1>  -  l  (n-lHx^j  "  xj)»'(x1) 


:  In  L*. 


denote  the  distinct  wnccnaored  failure  times  bv  <  T,  <  ...  *  T^, 
whcce  k  la  the  number  of  ^  which  equal  one.  Let  denote  the  number  of 
losses  which  occut  in  the  Interval  T^»  Including  any  loaaes  at 

but  not  at  T  J,  J  -  o.  1 . k.  where  TQ  -  0  and  -  -.  Let  the 

times  of  the  A  loaeea  be  denoted  bv  L^  ,  1  -  1»  *•  •••*  'j • 

Since  r(t)  la  increasing,  r(L^)  *  r(Tj),  1  “  1 . ^  (for  «*ch 

Xj  that  la  not  rero).  Therefore,  for  each  J,  0  S  J  <  k,  we  have 


V— V1 

-  I  ,  (n-l) (x,+j  -  »,)!(«,) 

1-X  *..M 


-  -in  -  (X0+...4l  1*J)}(L1l,)  -  T) 


)r  (T . ) 


-in  -  (»0*...u  -  t^htaj1’) 


-in  -  <»„♦.  -  tJjh.CLjj’) 
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-in  -  <V*‘**Vl+J)1V  r(V‘ 


replacing  -r(L^)  by  mo,  1*1 . 'Q-  Hence 


k  ^  ft* 

In  l*  s  I  In  r (T  )  -  I  a  r(T  )  -  In  L  , 
j-1  J  J-l  ’  1 


(2.4) 


where 


-in  -  (A 


j-e})  }Tj .  J  -  1.  •  •••  k~lt 


Now,  the  problea  ot  obtaining  the  naxlnuai  lilellhood  esclattor  of  r(t) 
aubject  to  the  condition  that  r(t>  1*  increasing  is  reduced  to  that  of 
max  Usl ring  In  L**  given  by  (2.4)  subject  to  the  constraint 

r(Tj)  s  r(T2)  s  ...  *  rty.  Let  v  -  rd^).  j  -  1 . V-  Then  we  with  to 


obtain 


k  k 

■ax  (  [  In  y  -  lav) 
J-1  3  J-1  3  3 

aubject  to  Tj  S  y,  S  ...  i 


(2.S) 


note  that  if  for  ao»e  J,  a^  -  0,  then  the  function  C(y^ . *1^ 


l  (In  y  -  a  y  )  la  not  bounded.  However,  when  J  4  k, 
J-1  J  J  J 
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a.  >  l,(3)-*-...+L.(J)*  ln-(*  ♦...♦A  .♦3))T. 

J  *  Aj  U  J  J 


-  +  .♦L,°)  -  *  T  5  0. 

J  3  3 


Therefore,  only  a.  can  be  zero  and  thia  occur*  when  there  are  no  censored 
obaarva 1 1  on*  larger  than  T^,  the  largest  uncensored  failure  tiAe.  If  a^  *  0, 
It  la  Impossible  to  obtain  a  maximum  likelihood  estia>ator  of  r(t)  directly 
by  solving  (2.3).  Consequently ,  we  first  consider  the  subclass  of 
distribution*  F  with  corresponding  (allure  rate  functlona  bounded  by  a 
constant  M.  Following  the  argument  of  Marstrall  t>  Prose han  (1963)  and 
utilizing  the  result*  of  Barlow,  Bartholomew.  Bremner,  4  Brunk  (197?,  p.  *4), 

Jal 

the  naxlaum  likelihood  estimator  of  r(t)  for  F  •  F  ‘  1*  given  by 


r^(T.)  •  min  [  min  max  <(v-u)/(r  3+...*r  3,)},  Nj 
n  1  vfHl  usi  U  V_1 


(2.6) 


where  r,  ■  M  and  r  •  ,  J  *  1,  2,  ....  k-1.  letting  M  -*  *’  In  (2.6),  we 

*  J  J 

obtain  the  maxinus  likelihood  estimator  of  r(t)  for  F  as 


where 


t  <  T. 


rn(Ti>.  T1  ?  *  Tl^r  1  "  l*  2 . k_1  (2,7) 


rn<T)l) ,  Tk  5  1 


r  (T. )  -  min  max  {(v-uKir  '♦...♦r  ,)).  1*1,  ....  k-1, 

°  1  vMM  usl  u  V'3 


and  r  (T.  )  "  •. 

n  r 


If  a^  i  0,  the  solution  of  (2.5)  Is  obtained  by  applying  the  results 
of  barlow,  Bartholomew,  Brenner,  4  Br  ink.  (1972,  p.  *4).  In  this  esse,  the 
maximum  likelihood  tst  lautor  of  r(t)  is  given  by  (2.7)  where 

*  -1  -1  ,  - 
r  (T  )  *  min  max  {(v-u)/(ru  + . . .♦ r y_ j ) } , 

°  v? 1+1  u ■.  1 

and  r^  •  a^*,  J  •  1,  2 . k. 

In  either  c*te,  the  maximum  likelihood  estimator  of  F(t)  is  obtained 
from  equation  (2.2)  as 

F  ( t )  ■  exp  {-/*  r(u)du) 
o 

•  exp  -I  r^(T(){mln  (t,  T^,)  -  T^)],  t  i  0,  (2.8) 

ii:  i?0,T  • t} 

where  •  0  and  •  m.  We  note  that  this  estimator  of  the  survival 

function  is  well-defined  lor  all  t  *•  0,  is  a  smooth  function,  and  approaches 
zero  as  t  ♦  »  (If  a.  •  0,  F(t)  •  0  for  t  i  T  since  r  (T  )  -  •>) . 

m  k  n  ® 

Similar  techniques  can  be  applied  for  the  case  that  F  has  decreasing 
failure  rate. 

3.  AN  EXAMPLE 

«'•  use  the  data  given  by  Kaplan  k  Meier  (1958,  p.  464)  and  also  used 
by  Susarla  4  Van  Byxln  (1976)  to  obtain  an  estimate  of  the  survival  function 
from  the  maximum  likelihood  procedure  In  Section  2.  The  ordered  data  are 
0.8.  1.0*,  2.7*,  3.1,  5.4,  7.0*,  9.2,  12.1*  months,  where  a  denotes  s  time 
of  lose.  In  our  notation  6  •  1,  1  •  1,  4,  5,  7  and  •  0,  1  ■  2,  3,  6,  8, 


with  T.  -  0.8,  T,  -  3.1,  T,  -  5.4,  and  T,  -  9.2.  Also,  a  -  0,  *  -  2, 

1  Z  j  *•  U  1 
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Lj4)  -  12.1.  The  a^  are  then  •  13.6,  a^  •  a}  •  9.2,  and  a ^  -  2.9. 

(Figures  1  and  2  about  here) 

Figure  1  shows  the  estimate  r0(t)  of  the  failure  rate  function.  We 
note  that  since  a  censored  value  was  observed  larger  than  T^,  r^d)  1“ 
finite  lor  all  t.  Figure  2  gives  the  estimates  of  the  survival  curve  using 
F(t)  given  by  (2.8)  and  Pvt’',  the  product  limit  estimate.  These  results 
can  be  compared  with  the  nonparanetrlc  Bayes  estimate  for  the  same  data 
given  by  Susarla  &  Van  Ryiin  (1976,  p.  900).  We  renark  that  to  obtain  the 
nonpat ametr lc  Bayes  estimate,  the  parameter  J  of  the  Dlrlchlet  process  prior 
must  be  chosen.  Susarla  &  Van  Ryzin  (1976)  indicate  the  effects  of  three 
choices  of  3  on  their  estimate. 

A.  SMAM  SA.HTLF  COMPARISONS 

We  have  performed  Monte  Carlo  simulations  for  the  Welbull  distribution 

in  order  to  obtain  the  email  sample  behavior  of  the  estimator  F(t)  as 

compsred  with  the  product  limit  estimator  P(t).  The  slmulstlons  were  based 

or  2000  random  samples  each  of  else  n,  { X ^ ,  from  a  Welbull  distribution 

-  a 

w* th  survival  function  F(t)  »  exp(-t  / 8 ) .  t  t  0,  with  right  cenaorshlp. 

The  censoring  random  variables  C, . U  were  chosen  to  he  Independent 

1  n 

of  the  x|  and  independent,  identically  distributed  as  uniform  on  (0,1^), 
where  T.  was  the  Cth  percentile  of  the  Welbull  distribution.  For  example, 
when  C  •  75,  we  obtain  on  the  average  251  censoring  in  the  samples  (censoring 
f  rsr  t Ion  0.25). 

The  average  squared  error  of  the  estimates  was  computed  from  the  2000 
trials  for  various  vslues  of  a  1,  n,  and  censoring  fractions.  Table  1  shows 


val  Probability. 


1  It 


(Table  1  about  here) 


loa«  of  the  result*.  Aa  anticipated,  the  performance  of  F  Improves  as  the 
censor  Ins  be. user  nore  t<vrir,  ai  a  Increases,  and  a*  n  Increases.  Compared 
with  P,  F  does  such  better  In  both  tails  of  the  distribution,  and  perform* 
a*  well  as  P  In  the  renter  based  on  the  mean  squared  error  .  This  Is  not 

surprising  for  the  upper  tall,  hovevei  ,  since  P  Is  not  we  1 1  -def  itied  when 

* 

the  largest  observation  Is  a  loss.  In  the  simulations,  we  defined  lit)  to 

be  the  constant  P(T^)  for  t  >  .  Ktsulis  similar  t»>  Table  1  were  also  found  when 

uncanaored  asm; let  were  used;  that  la.  the  estimator  of  Crenander  (1936) 

and  Marshall  fc  Prose  ha n  (1965)  behaves  In  the  same  manner  aa  Indicated  by 

Table  1  when  compared  with  the  product  limit  estimator  (which  Is  one  minus 

the  empirical  distribution  function  for  s  ccmpletc  sample). 
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Life  testing;  Tioduct  limit  estimator;  Right  censorship;  Nonpar are t ric 
estimation  of  failure  rate. 
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The  maximum  likelihood  estimator  of  a  distribution  function  with  monotone 
failure  rate  is  derived  based  on  a  sot  of  observations  subirct  to  arbitrate 
right  censorship.  This  estimator  is  defined  everywhere  on  the  positive 
real  line  while  the  Kaplan-Mcicr  estimator  rav  not  be.  The  small  sample 
properties  of  this  estimator  are  Indicated  by  results  of  a  Monte  Carlo 
studv  for  the  Wclbull  d  i  s 1 1 ibut ion . 
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obtained  for  triangular  arrays  of  rjndom  elements.  Finally,  the  direct 
applications  of  these  results  in  obtaining  consistency  of  the  kerney  den 
slty  estimates  are  Indicated. 
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